Asymptotic  Behavior  of  an  SI  Epidemic 
Model  with  Pulse  Removal 


K.M.  Fuhrman,  I.G.  Lauko  and  G.A.  Pinter 

University  of  Wisconsin- Milwaukee,  Department  of  Mathematical  Sciences, 

Milwaukee,  WI  53201-0413 


Abstract 

In  this  paper  we  discuss  an  SI  epidemic  model  with  pulse  removal  from  the  infective 
class  at  fixed  time  intervals  with  both  exponential  and  logistic  type  underlying  pop¬ 
ulation  dynamics.  This  model  has  a  significance  when  dealing  with  animal  diseases 
with  no  recovery,  or  when  we  consider  isolation  in  human  diseases.  We  provide  a 
rigorous  analysis  of  the  asymptotic  behavior  of  the  percentage  of  infected  individu¬ 
als,  the  total  number  of  infected  individuals,  and  the  total  population  in  our  model. 
We  show  that  periodic  removal/isolation  is  a  feasible  strategy  to  control  the  spread 
of  the  disease. 
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1  Introduction 


Mathematical  methods  are  widely  used  for  understanding  mechanisms  in  the 
spread  of  infectious  diseases.  Together  with  computer  simulations,  epidemic 
models  are  very  useful  in  explaining  and  testing  various  theories,  determining 
key  parameters,  and  finding  the  effects  of  changes  in  the  parameter  values. 
The  high  interest  can  be  explained  by  the  economic  and  societal  effects  of 
various  known  epidemic  diseases  and  the  emergence  of  new  viruses  and  in¬ 
fections.  Another  important  motivation  for  the  development  of  this  field  is 
the  evaluation  of  various  vaccination/control  strategies  for  human  as  well  as 
animal  diseases.  The  study  of  epidemics  among  animals  has  gained  a  lot  of 
attention  in  recent  years  with  major  outbreaks  of  the  mad  cow  and  hoof-and- 
mouth  diseases  together  with  the  spread  of  chronic  wasting  disease  in  the  US 

[9,10]- 

Infectious  animal  diseases  have  a  great  effect  in  the  ecology  and  management 
of  wildlife  and  domestic  animals.  Sometimes  the  diseases  among  animals  pose 
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a  great  threat  to  public  health.  When  there  is  an  epidemic  outbreak  among 
animals,  a  common  way  to  control  it  is  by  destroying  a  fraction  of  the  infected 
population.  In  [10]  Heesterbeek  and  Roberts  considered  a  model  with  a  fixed 
culling  rate  7  and  showed  its  effect  on  the  dynamics  of  the  infected  population. 
However,  the  constant  removal  strategy  assumes  continuous  culling  at  the 
same  rate,  which  is  not  a  realistic  approach.  We  propose  a  new  model  that 
assumes  that  a  fraction  of  the  infected  population  is  being  removed  at  certain 
fixed  points  in  time  and  not  continuously. 

The  effects  of  pulse  removal  have  been  studied  before  in  relation  to  epidemic 
models  while  analyzing  vaccination  methods  [1,2, 7, 8],  where  a  fixed  portion  of 
the  susceptible  population  was  removed  (i.e.  vaccinated)  at  fixed  time  inter¬ 
vals.  It  was  shown  that  the  infection  free  solution  can  be  obtained  with  lower 
levels  of  vaccination  compared  to  the  constant  vaccination  strategy.  Our  goal 
is  to  show  that  pulse  removal  from  the  infected  class,  which  can  model  the 
actual  destruction  of  a  portion  of  the  population  or  isolation  from  the  rest  of 
the  population,  is  a  feasible  strategy  to  control  the  disease.  The  control  ob¬ 
jective  is  to  drive  the  number  of  infected  individuals  to  zero  while  stabilizing 
the  population  size  or  growth  rate. 

In  section  2  we  consider  an  SI  model  with  exponential  underlying  dynamics 
for  the  total  population.  The  inclusion  of  pulse  removal  transforms  the  system 
into  a  system  of  impulsive  differential  equations  [4,5].  We  establish  a  threshold 
value  pi  for  the  removed  fraction  p  such  that  for  p  <  p\  the  fraction  of  infected 
individuals  tends  to  a  positive  periodic  function,  while  for  p  >  p\  this  fraction 
tends  to  0.  However,  the  fact  that  the  fraction  of  infected  individuals  tends  to 
0  does  not  guarantee  the  eradication  of  the  disease  in  a  growing  population. 
Thus,  it  is  of  interest  that  we  establish  another  threshold  value,  p%,  which  has 
the  property  that  for  p  >  p*2  the  number  of  infected  individuals  tends  to  0, 
while  the  population  grows.  This  case  realizes  the  control  objective.  We  also 
provide  a  detailed  analysis  of  the  asymptotic  behavior  of  the  total  population 
for  p  <  P2,  and  illustrate  the  different  possible  cases. 

In  section  3  we  carry  out  an  analogous  study  for  an  SI  model  with  logistic  dy¬ 
namics  for  the  total  population,  which  may  provide  a  more  realistic  approach 
[3,6,11,12], 


2  SI  Model  with  Pulse  Removal 


We  consider  a  common  SI  model  with  exponential  dynamics  for  the  total 
population.  When  pulse  removal  is  introduced  in  this  model,  we  obtain  the 
following  system: 
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(1) 

(2) 

(3) 


dS  JO  PSI 

dt  N  ’ 

-  =  -(d  +  a)I  +  ^  -  p  £  I(nT~)S(t  -  nT), 
S  +  I  =  N. 


Here  /(f)  >  0,  S(t)  >  0  denote  the  size  of  the  infected  and  susceptible  popu¬ 
lation  classes,  respectively,  and  N(t)  is  the  total  population  at  time  t  >  0.  b  is 
the  birth  rate,  d  is  the  death  rate,  while  (3  denotes  the  contact  rate  and  a  is  an 
added  death  rate  due  to  the  disease.  We  consider  the  case  when  the  natural 
birth  rate  b  is  greater  then  the  death  rate  d.  In  a  population  that  is  dying 
out  naturally,  that  is,  if  d  >  b,  the  removal  does  not  affect  the  asymptotic 
behavior  of  the  system.  Note  that  (l)-(3)  implies  that 


dN  °° 

——  —  (b  —  d)N  —  al  —  I{nT~)5(t  —  nT).  (4) 

dt  n= 1 

In  our  proposed  model,  the  removal  is  applied  as  an  impulse  whenever  t  =  nT, 
n  —  1,  2, ...  The  fraction  p  of  the  infective  individuals  is  being  removed  from 
the  class,  that  is,  /(f)  decreases  with  proportion  p  whenever  t  —  nT.  We  will 
use  the  notation 


I(nT  )  =  lim  /(f)  and  /(nT+)  =  lim  /(f)  (5) 

t—tnT-  t^nT+ 

to  indicate  the  different  limits  of  the  solution  from  left  and  right  at  these 
points.  The  sequence  of  impulses  creates  jump  discontinuities  in  /(f),  which 
decreases  with  proportion  p  whenever  f  =  nT.  We  start  the  analysis  by  trans¬ 
lating  the  model  for  infectious  and  susceptible  fractions  of  the  population. 


2. 1  Susceptible  and  Infective  Individuals  as  Fractions  of  Population 


We  translate  the  given  system  into  a  model  describing  the  behavior  of  s(f)  = 
§ffy,i(t)  =  Now,  s(t)  and  i(t)  represent  the  percentage  or  fraction  of  the 
population  that  are  susceptible  and  infective,  respectively. 


Thus  the  system  (l)-(3)  becomes: 


ds 

dt 

di 

dt 


b  -  bs  -  (P  -  a)si  +  f)  p'lfL  S(t  -  nT) 


n— 1 


1  —  pi{nT~ ) 


.  2  i(nT  )(1  —  i(nT  ))  _. 

—  (b  +  a)i  +  P  st  +  m2  —  ^2p - - - ~7IFFT — ~d{t  —  nT) 


n— 1 


1  —  pi(nT  ) 


(6) 

(7) 
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with  feasibility  set  V  —  {(s,  i)  \  s  >  0,  i  >  0,  s  +  i  —  1},  and  we  define 
T>i  —  V  —  {(1,0)}.  We  note  that  the  system  is  well-posed  epidemiologically, 
that  is,  if  the  initial  conditions  i0  and  s0  are  in  the  feasibility  region,  then  s(t ) 
and  i(t )  will  stay  in  V  for  all  t  >  0. 


Theorem  1:  Consider  the  system  of  equations  (6) -(7)  with  b  >  0,  f3  —  (b  + 
a)  >  0  and  all  other  parameters  non-negative.  The  following  holds  with  pi  = 

1  _  e-(p-b-a)T  . 

(1)  Up  >  p\  then  the  disease  free  equilibrium  ( s ,  i)  —  (1,  0)  always  exists  and 
it  is  globally  asymptotically  stable  in  the  feasibility  set  V. 

(2)  Ifp  <  pi  then  there  exists  a  unique  endemic  periodic  solution  ( i *,  s *)  and 
this  solution  is  globally  asymptotically  stable  in  V,. 

Proof:  Since  s{t)  —  1  —  i(t),  it  is  sufficient  to  consider  the  initial  value  problem 


di 

dt 


—  {b  +  a)i  +  (3(1  —  i)i  +  ai2 


~Ep 


i(nT~)(  1  -  i(nT~ )) 


n— 1 

i(0)  —io>  0. 


1  —  pi(nT~ ) 


S(t  -  nT) 


(8) 

(9) 


On  an  interval  [nT,  (n  +  1  )T]  (8)  is  a  logistic  equation  in  i  and  it  can  be 
explicitly  solved.  We  note  that  the  assumption  (3  >  (b  +  a)  guarantees  that  the 
equation  (8)  without  the  pulse  removal  has  two  stationary  solutions,  ii(t)  =  0 
and  i2(t )  =  /3~^a) ,  with  0  <  i2(t)  <  1-  The  trivial  solution  ii(t)  is  unstable, 
while  i2(t )  is  asymptotically  stable.  If  (3  <  (b  +  o),  then  there  are  two  cases: 


a)  If  (3  —  a  >  0,  then  i2(t )  <  0,  and  in  this  case  ii(t)  =  0  is  asymptotically 
stable,  in  particular,  for  any  i(0)  =  io,  0  <  io  <  1,  i(t )  — )-  0  as  t  — *  oo, 
where  i  is  the  solution  of  (8)  without  the  pulse  removal.  Thus  the  frac¬ 
tion  of  infected  individuals  goes  to  zero  naturally  without  any  controlling 
mechanism. 

b)  If  (3  —  a  <  0,  then  i2(t )  >  1,  and  we  again  have  that  for  any  0  <  io  <  1, 
i(t )  — *  0  as  t  —>  oo. 

We  also  remark  that  the  assumption  (3  >  (b  +  a)  is  realistic  in  epidemiological 
applications.  Solving  (8)  on  an  interval  [nT,  (n  +  1)T]  we  obtain 


G(t ) 


i(nT+)e^-b-a^t~nT^ 

l  I  i(nT+)(P-a)  ( Q(p -b-aMt-nT )  _  2) 
3—b—ot  ^  ' 


At  the  removal  points  we  have 
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i((n  +  1)T+) 


(l-p)G((n  +  l)r~) 

1  —  pG((n  +  1)T_)  ' 


Under  our  assumptions  this  solution  is  well-defined  and  exists  for  all  t  >  0. 


Now,  we  consider  the  sequence  of  the  “initial”  values  for  each  interval  between 
removals  at  i(nT+),  where  n  =  1,2,....  We  let  i{nT+ )  =  in  and  deduce  a 
stroboscopic  map  F  that  satisfies 


in+i  =  F{in).  (10) 

The  map  F  determines  the  value  of  i(t)  immediately  after  each  removal  at 
discrete  times  t  =  nT.  Then, 


*n+l  F{ir j) 


(1  -p)i{(n  +  1)T  ) 

1  -  pi((n  +  1)T~) 

(1  —  p)ine^~h~°^T 

1  +  -  1)  -  pineW-b~a)T' 


If  P  ^ 


{P-a)e^-b~a'>T-l 
(P—b—a)e^~b~a^T  ’ 


the  map  F  has  two  fixed  points 


*i  =  0  and  i2  = 


(1  —  p)e^  b  a)T  —  1 


Jt^_(e(P-b-a)T  _  !)  —  pe(P-b-a)T 


(11) 


We  analyze  i2  as  a  function  of  p.  Figure  1  illustrates  the  general  behavior  of 


g(p) 

1 

1 

1 

1 

1  _ 

(p-b-a)  - 

V 

(P-a) 

p 

(1_e-(P-b-a)T0)  | 

Fig.  1. 

General  behavior  of  the  function  g{p ) 

the  function  i2  =  g{p).  For  p  =  0,  we  have  g(0)  =  ^~^~a .  The  function  has 
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an  asymptote  at  pa  =  ’  w^ich  nray  or  may  not  be  less  than 

1.  The  function  is  always  decreasing  for  p  <  pa ,  and,  limp_^-  g{p )  =  —  oo. 
The  zero  of  the  function  is  at  p  =  p\  =  1  —  e~^~b~°^T  <  pa.  If  p  >  p\  then 
g(j>)  <  0  or  g{p)  >  1,  which  brings  it  outside  the  feasibility  region  D.  Finding 
the  derivative  of  F  at  the  two  fixed  points  will  determine  the  stability  of  those 
fixed  points.  We  obtain 

dF(in)  e(/3-6-a)T(f  _  p) 

din  =  (1  +  -  1)  -  pine^-»-)Ty ' 

The  derivative  at  the  fixed  points  yields 


dF(in) 

din 

dF{in ) 


(1  -p)e^-b~a)T, 

1 

(1  —  p'je^~b~a^T 


In  order  for  i\  —  0  to  be  locally  asymptotically  stable  we  make 


dF(in) 


di 


<  1 


which  implies 


p  >  1  -  e-(0-&-°)r  (12) 

If  condition  (12)  holds,  then  the  other  fixed  point  i*2  is  not  in  the  feasibility 
set,  that  is,  either  i2  <  0  of  i2  >  1.  The  above  argument  also  shows  that  in 
the  case  when  p  <  1  —  e_(^_6_a)T,  i2  is  locally  asymptotically  stable,  while  i\ 
is  unstable. 

Now,  we  show  the  global  stability  of  the  fixed  points.  Since, 


dF{in)  = _ b  a)T(l  -p) _ 

din  (1  +  ^in(e^-b~^T  -  1)  -  pine^-b-«)Ty 

for  any  0  <  in  <  1,  we  have  that  {Fn(i0)}  =  {i0,  F(i0),  F(F(i0)), ...}  is  a 
monotone  sequence  for  any  initial  condition  i0.  This  sequence  is  bounded  be¬ 
tween  0  and  1  by  the  well-posedness  of  our  system.  Therefore,  we  have  ob¬ 
tained  a  monotone  bounded  sequence  which  has  to  converge  to  some  limit  L. 
We  proceed  by  finding  this  limit. 
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L  =  lim  Fn(i0) 


(1  -p)Le 


(, 8-b-a)T 


j3—b—a 

Solving  for  L  we  obtain  the  following  result 


1  +  wE^L^~b~a)T  -  !)  -  pLe(P-b~a)T ' 


T  .-nr  f*  “  P)e{P  b  a)T  ~ 1 

!  “  h  “  ’  2  2  “  Jt_Z_(e(P-b-a)T  _  1}  _  pe(P-b-a)T  ' 

If  condition  (12)  holds  then  L2  ^  £>,,  so  =  0  is  the  limit  of  {Fn(i0)}  for 
any  i0,  and  thus  it  is  a  globally  asymptotically  stable  fixed  point  of  F(in). 
Hence  in  this  case  i(t)  — »  0  as  t  — >  oo  for  any  initial  condition  i0  e  2X  If 
p  <  1  —  e~^~b~a^T ,  then  L2  is  the  limit  of  (F”(f0)}  and  thus  L2  is  globally 
asymptotically  stable  in  V,.  Then  there  is  a  unique  endemic  periodic  solution 
iffi)  =  i(t )  with  i0  —  L2.  This  periodic  endemic  solution  is  globally  asymp¬ 
totically  stable,  meaning  that  |  i(t)  —  i^it)  |— >•  0  as  t  — »  oo  for  any  i(t)  which 
solves  (l)-(3).  This  concludes  the  proof  of  Theorem  1. 


2.2  Total  Number  of  Infective  Individuals 


Now,  we  generalize  the  result  from  the  previous  section  to  study  the  behavior 
of  the  total  number  of  infective  individuals  and  the  total  population.  We  go 
back  to  the  system  of  equations  (l)-(3). 

Theorem  2:  Consider  the  system  of  equations  (l)-(3)  under  the  assumption 
that  b  >  d  >  0,  (3  —  (ck  +  6)  >  0  and  all  other  parameters  non-negative.  Then 
with  P2  —  1  —  e_^_d_a)T  the  following  hold: 

(1)  If  p  >  1-  e-C-d-o)T  the  total  infected  population  has  the  asymptotic 
behavior  I(t)  0. 

(2)  If  1  —  <  p  <  1  —  e~^~d~a')T  the  total  infected  population  has 

the  asymptotic  behavior  I(t)  — »  00. 

Proof:  First,  we  consider  the  case  when  p  >  1  —  e~^~b~a^T  ^  that  is,  the 
fraction  of  infected  population  goes  to  0.  By  the  previous  theorem,  we  know 
that  under  this  condition,  s(t)  =  — y  1.  Then  for  any  s  >  0  there  exists  t0 

such  that  1  —  s  <  jiM-  <  1  for  t  >  to.  Let  to  =  (noT)+.  Consider, 


dT  BST  00 

—  =  -(d  +  a)I  +  f-j^--pJ2l(nT-)6(t-t0-nT),  t  >  t0  (13) 

with  I  (to)  =  Jo-  Under  the  assumption  that  1  —  £  <  <  1  for  t  >  to,  we 

introduce 
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(14) 


dl  -  00  „ 

—  =  —  (d  +  a)I  +  (3(1  —  e)I  —  p  ^  I(nT~)5(t  —  to  —  nT) 

dt  n=  1 

with  I  (to)  —  Io- 

Then,  I(t )  >  I(t )  for  all  t  >  to  and 

/(f)  =  /0[(1  _  p)em-s)-d-<*)T}nem-E)-d-a)(t-nT)  ^ 

for  nT  <  t  <  (n  +  1)T,  n  >  n 0. 

Thus,  the  sequence  of  ’’initial”  values  at  I(nT+ )  where  n  =  no,  no  + 1,  no  +  2... 
forms  a  geometric  progression,  which  converges  to  0  if  (1  —  p)e(/3(1~£'>~d~a')T  <  1 
and  grows  without  bounds  in  the  reverse  case. 

Therefore,  we  conclude  that  if  p  <  1  -  e-(£(l -e)-d-a)T  then  I(t )  — >  00  as 
f  — »  00.  Since  /(t)  >  /(f)  for  f  ^  to  and  s  0  is  arbitrary,  W6  have  I (t^j  — y  oc 
as  t  — >  oo  if  p  <  1  —  e_(^_d_a)T . 

Next  let  us  consider 


dJ  00  _ 

—  =  -(d  +  a)/  +  (31  -  p  Y,  I(nT~)5(t  -  f0  -  nT) 

dd  n= 1 

with  I (t0)  —  lo¬ 
in  this  case,  we  have  I(t)  <  /(f).  By  similar  analysis  we  obtain  that  for 
p  >  1  —  e~^~d~a',T  we  have  /(f)  — >  0  as  f  — >  oo.  Since  0  <  /(f)  <  7(f),  we 
have  that  I(t)  — f  0  as  well,  which  concludes  the  proof  of  Theorem  (6.2). 

We  note  that  the  assumption  b  >  d  guarantees  that  the  threshold  pi  for  /(f) 
is  greater  than  p\. 

Now,  we  consider  the  behavior  of  /(f)  when  p  <  1  —  e~^~b~°l')T .  In  this  case, 
i(t)  approaches  an  endemic  periodic  solution  il(t),  for  which  i^(nT+)  =  i*2, 
where  i*2  is  the  nonzero  fixed  point  of  the  map  F  in  (10)  .  We  introduce 


H(p)  =  (1  -  p)e^-d-a)T(/  1  a  i*2(e^-b~a)T  -  1)  +  1)^  (15) 

(3  —  b  —  a 

where,  i*2  is  given  by  (11).  The  definition  of  H(t)  is  motivated  later  in  the 
proof  of  the  following  theorem,  which  describes  the  dynamics  of  /(f)  under 
the  conditions  p  <  1  —  e~^~b~a^T  and  (3  >  (b  +  a). 
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Theorem  3:  Consider  the  system  of  equations  (l)-(3)  with  b  >  d  >  0,  /3  > 
(i b  +  a)  and  all  other  parameters  non-negative.  Let  p  <  1  —  e~^~b~°^T .  In  this 
case  we  have: 

(1)  If  H(p)  <1,  then  the  total  inf ected  population  has  the  asymptotic  behavior 

/(f)  ->  0. 

(2)  If  H(p)  >  l,  the  total  infected  population  has  the  asymptotic  behavior 
I(t)  — >  oc. 

Proof:  From  (2)  we  have  that: 


dJ  B3T  °° 

-  =  -(i+c)I+P—-p  ■£nnT-)5(t  -  nT). 

Substituting  sit)  —  S(t)/N(t)  we  obtain 


dl  °° 

—  =  (3s(t)I  -(a  +  d)I-p^2  I(nT~)S(t  -  nT ), 

dt  n= 1 


then 


dl 

dt 


=  (13(1  -  i(t ))  -  (a  +  d))I  ~pJ2  I(nT~)S(t  -  nT). 


n— 1 


We  solve  this  equation  for  (n  —  1)T  <  t  <  nT  and  obtain 


I(t)  =  I0efo 


/3(1 —i(r))—(a-\-d)dT 


,  where  Iq  —  I((n  —  1)T+). 


From  Theorem  1,  we  know  that  i(t)  approaches  an  endemic  periodic  solution 
i^(t),  with  i*2(nT+)  —  i*2.  Thus,  i(r)  is  approximately 


(n  —  1  )T  <  t  <  nT  for  large  n. 


i+«&*2(e(/3“6“a)T-1): 


We  consider  I(nT+),  the  number  of  infective  individuals  right  after  the  re¬ 
movals.  Then  we  obtain  the  following: 


I(nT+)  =  (l-p)I0e 


09 -d-a)T 


fT  0 

Jo  -P- - 


0-b-a 2 


,;e(d-l-a)r+(1.?^2.j.) 


dr 


=  (1  °  -  1)  +  l)s4. 

p  —  b  —  a 


Hence,  the  function 
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H(p)  =  (1  ~  p)e^~d~a)T ( ^  -  a  Jl(e^-b~a)T  -  1)  +  1)f£ 

determines  the  behavior  of  /(f).  The  above  calculation  shows  that  I(nT+)  = 
IoH(p),  where  Jo  =  I((n  —  1)T+).  Repeated  application  of  this  relationship 
yields 


I(mT+ )  =  I0(H(p)Ym-^-^  m>  n. 

Hence,  if  H(jp)  1  for  a  particular  p,  p  <C  1  —  e  ■■  ^  ^  a^,  then  /(f)  — y  oo  as 
f  — »  oo  for  this  value  of  p;  and  if  H(p)  <  1  then  /(f)  — *  0. 

In  the  next  section  we  show  the  examples  for  different  possible  behaviors  of 
this  function. 


2.3  Total  Population 


We  presented  a  complete  analysis  of  the  behavior  of  the  fraction  of  infective 
individuals  i(t)  and  the  total  infective  population  /(f).  Using  this  information 
we  can  determine  the  dynamics  of  the  total  population  N(t)  for  various  values 
of  p.  The  dynamics  of  the  solutions  of  (l)-(3)  can  be  summarized  as  follows 
for  various  values  of  the  parameter  p. 


p 

i(t) 

I(t) 

N(t) 

p  <  pi  and  H(p)  <  1 

i  — *  endemic  periodic  solution 

/  ->■  0 

N  ->■  0 

p  <  pi  and  H{p)  >  1 

i  — *  endemic  periodic  solution 

/  — >  oo 

N  — >  oo 

P*1<P<P*2 

i  0 

/  — )•  oo 

N  — >  oo 

P>P*2 

i  0 

/  ->■  0 

N  — f  oo 

Table  1:  Dynamics  of  solutions  of  (l)-(3)  for  various  values  of  parameter  p, 

with  b  >  d,  (3  >  (b  +  a) 

Therefore,  we  have  obtained  a  complete  description  of  the  asymptotic  behavior 
of  the  solutions  to  our  proposed  model. 

As  noted  before,  based  on  the  behavior  of  the  function  H(p)  we  can  have 
different  scenarios.  The  following  examples  depict  some  interesting  cases.  We 
note  that  the  parameter  values  in  the  examples  are  for  illustration  purposes, 
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they  do  not  necessarily  represent  real-life  situations.  We  use  the  following 
notation  for  the  threshold  values  of  p  in  the  interval  [0, 1] 

Phx  —  value  of  p  where  H(p)  —  1  and  H(p)  is  decreasing, 

Ph2  —  value  of  p  where  H(p )  —  1  and  H(jp)  is  increasing, 
p\  —  1  —  e~^~b~a')T ,  the  threshold  of  the  fraction  of  infected  population  z(f), 
p*2  —  1  —  e_^_d_a^T,  the  threshold  of  the  total  infected  population,  /(f) 


Example  1:  We  consider  the  case  when  first  H(p)  <  1  and  then  H(p)  >  1, 
hence  pi n  does  not  exist  in  [0,1].  This  scenario  can  be  described,  for  exam¬ 
ple,  by  the  following  parameter  values:  (3  =  2.25,  b  =  0.4,  d  =  0.1,  a  = 
0.8,  and  T  =  1.  The  thresholds  for  this  example  are  as  follows:  p/ 12  — 
.6197,  pi  —  .650062,  and  p*2  —  .74076.  First,  we  graph  the  case  when 
p  —  .5  <  Ph2.  Hence,  we  should  obtain  that  i  converges  to  the  endemic  periodic 
solution,  /  — >  0  and  N  — *  0.  The  results  of  the  numerical  simulation  for  this 
case  are  depicted  in  Figure  2. 


TIME 


Fig.  2.  Example  1  for  p  =  0.5 
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Thus,  we  can  see  that  the  larger  pulse  removal  p  in  this  case  has  the  effect  of 
’saving’  the  population  in  the  sense  that  now  N(t)  — *  oo  as  t  — *  oo.  At  the 
same  time  I(t)  — »  oo  also,  i.e.,  the  disease  is  sustained  in  the  population  as 
shown  in  Figure  3. 

If  ith  <  P  —  .7  <  pi,  then  i  — >  0,  I  — >  00,  and  A/-  — >  00.  The  larger  percentage 


Fig.  4.  Example  1  for  p  =  0.7 

of  removed  infective  individuals  now  results  in  the  fraction  of  the  infected 
population  i(t)  — >■  0,  while  J(f),lV(f )  — »  00  (see  Figure  4). 

If  p  —  0.77  >  P2  then  i  — >  0,  I  — >■  0,  and  N  — >  00.  Figure  5  shows  that  the 
removal  of  77%  of  infected  individuals  drives  the  total  number  of  infective  in¬ 
dividuals  to  zero,  while  the  population  is  growing  without  bounds.  Thus,  with 
this  control  strategy  the  disease  is  eradicated  while  population  is  flourishing. 
(This  is  in  contrast  with  the  no-control  (p  =  0)  case  where  the  population  is 
eradicated  together  with  the  disease.) 


3  SI  model  with  logistic  dynamics 


In  this  section  we  study  an  SI  model  with  pulse  removal  that  has  a  logistic 
underlying  dynamics  for  the  total  population  S  + 1.  This  is  considered  a  more 
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TIME  TIME 


Fig.  5.  Example  1  for  p  =  0.77 

realistic  approach  than  the  exponential  dynamics  especially  for  longer  time 
intervals.  We  consider  the  model 

^-  =  bN-d(l  +  rN)S-^-  (16) 

dt  N  y  ' 

RST  °° 

—  =  —  (d(  1  +  rN )  +  a) I  +  -  p  I(nT~)S{t  -  nT)  (17) 

N  —  S  +  I.  (18) 

This  implies  that  the  dynamics  of  the  total  population  is  given  as 

dN  00 

— —  =  {b  —  d)N  —  rdN2  -aI-pJ2  I(nT~)S(t  -  nT).  (19) 

dt  n= 1 

The  logistic  dynamics  is  the  result  of  the  assumption  that  the  death  rate  is  an 
increasing  function  of  the  total  population  in  the  form  d(  1  +  rN),  where  r  is 
a  small  positive  parameter.  The  following  theorem  describes  the  asymptotic 
behavior  of  the  infective  and  susceptible  classes. 

Theorem  4:  Consider  the  system  of  equations  (16)-(18)  with  b  >  0,  /3—  (b  + 
a)  >  0  and  all  other  parameters  non-negative.  We  have: 
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(1)  If  p  >  1  —  =  pi  then  the  disease  free  equilibrium  (s,i)  — 

(1,  0)  always  exists  and  it  is  globally  asymptotically  stable  in  the  feasibility 
region  V.  Moreover,  lit )  — >  0  and  N(t )  ^  as  t  — »  oo. 

(2)  Ifp  P\  then  there  exists  a  unique  endemic  periodic  solution  (  i  ,  s  )  and 
this  solution  is  globally  asymptotically  stable  in  V,.  In  this  case  I(t )  and 
N(t )  both  converge  to  a  positive  periodic  function,  so  the  disease  stays 
endemic  in  the  population. 

Proof:  We  proceed  as  in  Section  2  and  rewrite  equations  (16)-(18)  in  terms 
of  fractions  of  infected  and  susceptible  classes.  Thus  we  obtain 


ds  JS,  iinT  )s(nT  ) 

b(l  —  s)  —  {f3  —  a) si  +  p- “  nT ) 

n=  1  i 


—  pi{nT  ) 


dt 

u — x 

di  ,2  ^  i(nT~)il  -  i[nT~)) 

—  =  - (b  +  a )i  +  [3si  +  ai2  -  - ■(  - LS(t  -  nT) 

dt  n=i  1  —  piynT  ) 


with  feasibility  set  V  —  {(s,  i)  \  s  >  0,  i  >  0,  s  +  i  —  1.} 


(20) 

(21) 


We  can  see  that  equation  (21)  is  the  same  as  the  one  governing  the  behavior  of 
i(t)  with  exponential  underlying  population  dynamics  (7).  Hence,  our  previous 
results  are  directly  applicable.  Thus  we  have  that  for  p  >  1  -  e“^“6_a)  =  pi 
i(t )  — >  0  as  t  — >  oo,  while  i(t )  tends  to  a  positive  periodic  function  when 
p  <  p\.  Next  we  examine  the  behavior  of  I(t).  The  logistic  dynamics  of  the 
total  population  in  (19)  guarantees  the  boundedness  of  N,  so  i(t )  — >  0  implies 
that  the  total  number  of  infective  individuals  I(t)  tends  to  0  as  well.  This  in 
turn  implies  that  N(t)  — >  ^  (the  positive  equilibrium  of  (19)  with  /  =  0)  as 
t  — >  oo.  Thus  in  this  case  we  have  only  one  threshold  quantity  that  separates 
the  asymptotic  behavior  of  the  infected  and  susceptible  classes.  In  the  case 
when  p  <  p{  we  have  that  i(t)  tends  to  a  positive  endemic  solution  i*(t). 
Plugging  in  /(f)  =  i*(t)N(t)  into  (19)  yields  that  N(t)  also  tends  to  a  positive 
periodic  function.  Figures  6  and  7  illustrate  the  behavior  of  the  system  for 
p  >  pi  and  p  <  pi . 


Fig.  6.  Logistic  example  for  p  >  p\ ,  infected  class  (left),  total  population  (right) 
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Fig.  7.  Logistic  example  for  p  <  p\,  infected  class  (left),  total  population  (right) 
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